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ABSTRACT 

We present results of a numerical integration of the hydrodynamical equations govern- 
ing the self-similar, two-dimensional gas flow behind the bow shock of an FRII radio 
source embedded in an IGM with a power law density profile. The model predicts 
pressure gradients within the cocoons consistent with modest backflow. For very steep 
external density profiles sources may well not expand in a self-similar fashion and in 
this case the model is not self-consistent. The assumption of ram pressure confine- 
ment of the cocoons perpendicular to the jet axis is found to overestimate the ratio 
of the pressure in front of the radio hot spots and that in the cocoons. Based on the 
properties of the gas between bow shock and cocoon we calculate the X-ray surface 
brightness of the flow. This emission is found to be a good tracer of the density dis- 
tribution within the flow and varies significantly with the properties of the unshocked 
IGM. The cooling-time of the shocked IGM is found to be comparable to, or greater 
than, the Hubble time. The influence of a radio source on the evolution of its gaseous 
surroundings therefore extends well beyond the limited life time of the source itself. 
We compare our results with the X-ray map of Cygnus A and find some evidence for 
cold, dense gas clumps in the surroundings of this object. The extended X-ray emission 
observed around 3C 356 may also be caused by the bow shock of this radio source. 
We also present an empirical model for the X-ray emission of the shocked IGM due 
to thermal bremsstrahlung. 

Key words: hydrodynamics - galaxies: active - galaxies: individual: Cygnus A - 
galaxies: individual: 3C 356 - galaxies: jets - X-rays: galaxies 
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1 INTRODUCTION 

Most analytical models of FRII extragalactic radio sources 
are based on the existence of two jets emerging from an AGN 
in the core of the source; a model first proposed by Rees 
(1971). After passing through a strong shock at the location 
of the radio hot spot, the jet material inflates a cavity of hot, 
rarefied gas. The expansion of this cocoon is supersonic with 
respect to the surrounding IGM and therefore drives a bow 
shock into this material. This shock will compress and heat 
the gas of the IGM, thereby increasing X-ray emission from 
thermal bremsstrahlung. Inferred on theoretical grounds by 
Scheuer (1974), this hot layer of gas between bow shock 
and cocoon was first observed in Perseus A by Bohringer 
et al. (1993) and in Cygnus A by CarilU et al. (1994). The 
extended X-ray emission detected in the vicinity of 3C 356 
(Crawford & Fabian 1996) may represent another example of 
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thermal bremsstrahlung of the IGM boosted by the presence 
of a powerful radio source. 

The first analytical investigation of a cavity expanding 
supersonically into a surrounding atmosphere was presented 
by Sedov (1959) who considered the case of a strong explo- 
sion in a uniform density profile. Dyson et al. (1980) ex- 
tended this analysis to spherically symmetric winds with a 
permanent energy input and Heinz et al. (1998) included 
density gradients by modeling the external density distribu- 
tion with a King (1972) profile. The solutions of Sedov and 
Dyson et al. are both self-similar. The expansion of both 
the cocoons and bow shocks of FRII sources is shown to 
be self-similar regardless of its geometrical shape by Kaiser 
& Alexander (1997) (hereafter KA), and this result suggests 
that a solution based on the spherically symmetric case may 
be possible for a radio source provided that the elongation 
of the cocoon along the axis of the jet can be accommo- 
dated. The solution for the flow of the shocked IGM be- 
tween bow shock and cocoon requires solution of at least 
a two-dimensional problem instead of the one-dimensional 
spherical case of an explosion or wind. 

In this paper we assume the bow shock surrounding 
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the cocoon of FRII radio sources to be elliptical. We then 
determine the properties of the shocked gas between bow 
shock and cocoon by integrating numerically the relevant 
hydrodynamical equations. Analogous to KA, we will as- 
sume that the density distribution of the IGM surrounding 
the radio source is well described by a power law of the form 
Po{r/ao)~ , where r is the radial distance from the centre 
of the distribution. 

Although we may choose the shape of the bow shock 
freely without violating the conditions for self-similar expan- 
sion of the cocoon or bow shock, we will show that another 
condition for self-similarity may not be fulfilled by some of 
the solutions presented here. This condition is the pressure 
equilibrium of the jets within the uniform pressure of the 
cocoon away from the hot spots. This implies that for these 
cases the model may not be self-consistent. One possibility, 
as we discuss, is that the assumption of an elliptical shape 
for the bow shock must be relaxed and the source expands 
self-similar ly with a somewhat modified shape. However, we 
believe that the results obtained in this paper are an im- 
portant first step in solving the problems posed by the flow 
between bow shock and cocoon of FRII sources. 

The general dynamics of the large-scale structure of 
FRII sources depends on the details of the confinement of 
the cocoons of these objects. It is usually assumed that the 
pressure of the material in the cocoon is balanced by the ram 
pressure of the IGM (e.g. Scheuer 1974, KA). In this paper 
we show that the assumption of ram pressure confinement of 
the cocoons overestimates the ratio of the pressure in front 
of the hot spot and that of the cocoon material, leading to 
lower hot spot advance speeds in the analysis of KA. 

The results of the integration of the hydrodynamical 
equations allows us to determine the X-ray surface bright- 
ness of the shocked material. The result of this calculation 
can be directly compared with X-ray observations. The only 
currently available X-ray map of a powerful FRII source 
with sufficient spatial resolution to allow such a compari- 
son, is the ROSAT map of Cygnus A presented by Carilli et 
al. (1994). However, the advent of the new X-ray telescopes 
such as AXAF will allow us to extend this analysis to a 
greater number of FRII sources. 

In Section 2 we transform the hydrodynamical equa- 
tions to an elliptical coordinate system comoving with the 
self-similar expansion of the bow shock and the cocoon. Sec- 
tion 3 describes the numerical method used to integrate the 
equations. The properties of the flow are discussed in Section 
4. In Section 5 we calculate the X-ray emission as predicted 
by the model and compare our results with the X-ray maps 
of Cygnus A and 3C 356. 



2 THE HYDRODYNAMICAL EQUATIONS IN 
AN ELLIPTICAL COORDINATE SYSTEM 

The hydrodynamical equations governing adiabatic fluid 
flow in three dimensions if gravitational fields and viscos- 
ity can be neglected, are the equation of motion: 



— v-f (v-V)v = — Vp, 

Ot p 

the equation of continuity: 
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and the adiabatic equation: 
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where v is the velocity vector of the gas, p the pressure, 
p its density and F its adiabatic index. We assume in the 
following analysis that the bow shock surrounding the co- 
coon is a prolate ellipsoid with rotational symmetry about 
the 2-axis defined by the jet, and therefore we rewrite these 
equations in prolate spheroidal coordinates (ry, ^, <j) = tt/2, 
see Appendix |A.l.l[ ). The rotational symmetry reduces the 
problem to two spatial dimensions. 
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where all dashed quantities are functions of the ellip- 
tical coordinates rj and ^, c is constant and c,. = 
(sinh^ j?-|-sin^^)^''^. 

As the bow shock expands it is possible for all times 
to find a stationary prolate spheroidal coordinate system 
characterized by a constant c in which the bow shock is 
described as a coordinate surface with constant 77 = 775. The 
aspect ratio of the cocoon, i.e. the ratio of its length to its 
width, Ri,, and the 'length' of the bow shock, i.e. the point 
at which the cut of the bow shock with the a^z-plane crosses 
the z-axis, determine the particular set of values for c and 
rji,. KA show that the expansion of the bow shock and of 
the cocoon is self-similar; in this case the bow shock can be 
described as a coordinate surface with rj = rjb for all times. 
This is achieved by setting c — CoLt, where Co is a constant 
and Li, = f{t) is the 'length ' of the bow shock as described 
above (see Appendix A. 1.2). It follows that 
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Because of the self-similar evolution of both the bow 
shock and the cocoon the fiow between them must also be 
self-similar. This allows a further simplification of equations 
(g) by introducing dimensionless variables for the velocity, 
density and pressure of the shocked IGM in which the de- 
pendencies on time (i.e. on Lt and Lb) and the spatial co- 
ordinates are separated. 
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where we model the external density distribution with a 
power law, Po{r/ao)~ with po the central density, Uo the 
core radius and r the radial distance from the centre of the 
distribution. 

Using the comoving coordinate system described above 
and replacing the flow variables with their dimensionless 
counterparts we find 
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where we have split the equation of motion into two equa- 
tions corresponding to components parallel to e,, and e^ 
respectively, and 



O'q — ^ r] Cq Cj- Cg , Of — (vf 'T' Cq Cj- Cg . 
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Furthermore, the ratio of the length of the bow shock. 
Lb, and the length of the cocoon, Lc, must also be a constant. 
Using the result 



L, = ci{poa't) ^^ QT^ t^ 



(9) 



from KA, where Qo is the jet power and c\ is a dimension- 
less constant, LbLb/Ll reduces to (/9 — 2)/3, and the time 
dependence in equations (|7|) can be eliminated. 



3 NUMERICAL METHOD 

The set of equations (Fl) represent a system of first order 
partial differential equations. Although there are no time 
derivatives in this set the general form is that of an initial 
value problem because the values of the variables are known 
at the surface of the bow shock and the solution has to 
be propagated from this surface at rjb inwards to smaller r) 



in the direction of the contact discontinuity delineating the 
cocoon. 

The projection of the bow shock onto the a;z-plane is 
an ellipse with rib = constant and is described by 



Tb = 'i^\ Co sinh rjb + tanh rjb l\ 



(10) 



where lb is the dimensionless coordinate defined by lb — 
z/Lb- At the tip of the bow shock where h = 1, rb — 0, and 
for lb — 1/2 one finds Vb — 1/(2 Rb). These conditions yield 
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To derive the initial conditions at the bow shock we 
need to know the velocity of the bow shock perpendicular 
to its surface, uj_ , and from the appendix of KA we find 






. dz 



(12) 
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Since the bow shock expands self-similar ly, Dj, = xjLb and 
Vz = lb- By assuming that the bow shock is strong and that 
the adiabatic index of the external atmosphere is 5/3 this 
implies for the assumed elliptical shape of the shock surface 



R{vb,0 



3 sinh?;;, 

4 c, ' 
0, 






4Co (cos^ 


^ -1- sinh^ 


^ft) 


-Co*" (cos 


^ ^ -1- sinh" 


Vb) 



sinh^ rjb 



(13) 



Given the symmetry of the problem it is sufficient to in- 
tegrate the equations only between ^ = and ^ = 7r/2 which 
corresponds to the z and a;-axes respectively, and the follow- 
ing boundary conditions must apply along these coordinate 
lines: 



t/e = 0, 

|P - 0- (14) 

The numerical integration was performed using a finite 
difference discretisation scheme employing NAGP] routines 
on a grid of 1001 points between ^ = and ^ = n/2 sepa- 
rated by equal intervals of length 7r/2000. Equations (H) are 
non-linear and it is therefore not straightforward to assess 
the stability of the numerical method used. Trial calcula- 
tions show that in particular the solution for U^ is subject 
to oscillations. To suppress these oscillations, it is neces- 
sary to introduce an 'artificial viscosity' term of the form 
ad^U^/d^^ to all equations in the system containing ex- 
pressions proportional to dU^/d(, (e.g. Richtmyer & Morton 



T NAG is a registered trademark of Numerical Algorithms Group 
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1967). The constant a is chosen to be small compared to any 
other term in the system of equations in order to keep the 
influence of the 'artificial viscosity' within the error of the 
integration method. 

The integration has to be stopped at the contact dis- 
continuity. There can be no gas flow across this surface and 
therefore the criterion for stopping the integration is met 
when the component of the gas velocity perpendicular to 
the contact discontinuity is equal to the corresponding com- 
ponent of the self-similar expansion velocity of this surface. 
The shape of the contact surface, rc{z), and therefore the 
self-similar expansion velocity, ze^ + rc{z) e^, is not known 
a priori. Since the bow shock will be closest to the cocoon at 
the hot spot, one expects the fiow to first reach the expan- 
sion velocity of the cocoon at ^ = and some value rj = rjc, 
i.e. on the z-axis. Assuming that dvc/dz — > cx3 at this point, 
the condition for stopping the integration on the z-axis be- 
comes UnirjcS, — 0) — Co cosh 77c, where we have again made 
use of the fact that the cocoon is expanding self-similarly. 
A comparison with equation (H) shows that at this point 
6^ = which means that the system of equations no) breaks 
down, as one would expect at a contact discontinuity. 

The integration was advanced from the bow shock at 
r; = rib to smaller values of rj in steps of ry^/S x 10 and after 
each step it was checked whether the flow had reached the 
cocoon at the hot spot, i.e. whether the velocity of the flow 
along the z-axis was within 1% of the self-similar velocity. 

The cocoon is not elliptical and therefore does not co- 
incide with the coordinate surface at rjc- The integration 
was continued by considering the grid point of the finite dif- 
ference mesh of smallest ^ for which the solution had not 
yet reached the contact surface. All derivatives with respect 
to ^ in equations (|7|) were approximated by linear extrap- 
olation using two values of these expressions separated by 
775/5 X 10^ along a line of constant ^. This extrapolation effec- 
tively converts equations uo) into a system of ordinary differ- 
ential equations which were solved using an implementation 
of the backward differentiation formulae from the NAG ' li- 
brary. The integration was continued until the solution for 
the last point of the grid had reached the cocoon. 



4 THE GAS FLOW BETWEEN BOW SHOCK 
AND COCOON 

An important test of the analysis and integration presented 
here is to consider the case of a spherical symmetric wind as 
discussed by Dyson et al. (1980). These authors considered 
only the case of a wind expanding into an uniform atmo- 
sphere, but using the results of the previous section it is 
straightforward to extend their analysis to more complex 
environments in which /3 7^ (for a solution using King 
density profiles see Heinz et al. 1998). 

In the case of a spherical bow shock there is an addi- 
tional degree of rotational symmetry about the a;-axis. This 
implies that one can set ^ — 0; U^ and all derivatives with re- 
spect to ^ vanish. Furthermore it is possible to change from 
the remaining prolate spheroidal coordinate 77 to the more 
familiar spherical polar coordinate r by using r — Co cosh 77. 
Because of these simplifications the second equation in (H) 
vanishes completely and we obtain 
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where we have used Rb = l/\/3 for a spherical bow shock 
and therefore 77;, — > cxa which implies that coth 77 — > 1 for all 
77 close to 77;,. Setting j3 = recovers the equations used by 
Dyson et al. (1980). The initial conditions at the bow shock 
at 7- = 1 for the spherical case are U,- — 3/4, 7? = 4 and 
P = 3/4. 

Figure ^ shows plots of the solution for the case of a 
spherical bow shock. The curves for (3 — are identical 
to those presented by Dyson et al. (1980). In this case the 
contact discontinuity is reached at re = 0.86. For the two 
other cases presented here we find r^ = 0.85 and re = 0.84 
for 13 = 1 and /3 = 2 respectively. In all three cases the gas 
is accelerated away from the cocoon but the acceleration is 
somewhat less effective for the cases with a density gradient 
in the external atmosphere, hence the greater detachment of 
the cocoon from the bow shock. While this difference in the 
flow pattern is not very significant, the distribution of the 
density and the pressure of the gas behind the bow shock 
is distinctively different for each case. The gas is expanding 
adiabatically after being compressed by the bow shock. This 
behaviour is evident from the curves in Figure [ij for /3 = 
but is masked for /3 = 2 by the presence of dense gas at 
high pressure close to the cocoon which has passed through 
the bow shock earlier than the material in the regions closer 
to the bow shock. The density and pressure of the external 
gas in front of the bow shock decreases faster than that 
of the material which is expanding after passing through 
the bow shock, giving the impression that the gas is further 
compressed after passing through the bow shock. In this case 
the cocoon is surrounded by a relatively thin, heavy shell of 
highly compressed gas, while for (3 = the gas closest to the 
contact discontinuity is rarefied. 

Results from the numerical integrations of equations (r7|) 
for the more complicated case of ellipsoidal bow shocks ap- 
plicable to FRII sources are summarized in Table y. Note, 
that the stand-off distances between bow shock and cocoon 
are small both at the hot spot and close to the x-axis. The 
shape of the cocoon is therefore also close to an ellipsoid and 
this agrees well with observations (e.g. Leahy & Williams, 
1984). The following sections describe the results in greater 
detail. 



4.1 The velocity field 

Figure [4.l| shows a cut defined by the xz-p\&ne of the ve- 
locity field of the gas fiow between bow shock and contact 
discontinuity for one half of a wide {Rb = 1) ellipsoidal bow 
shock in an uniform atmosphere. Because of the rotational 
symmetry about the z-axis the vectors are lying fully in the 
xz-plane and there is a reflection symmetry about the z- 
axis. The vectors show the velocity of the flow with respect 
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to the comoving coordinate system, with the bow shock and 
cocoon at rest. 

A closer examination of the flow close to the hot spot 
shows the flow bending around the tip of the cocoon and 



then moving along the contact discontinuity (Figure 4.1 
left). On the z-axis, the axis of the jet, the flow ends in 
a stagnation point at the hot spot, similar to the case of a 
spherical bow shock. Close to the a;-axis the flow is decel- 
erated and there is a stagnation point on the x-axis ahead 
of the cocoon (Figure ]A^, right). This behaviour is imposed 
on the solution by the assumed symmetry. The flow pattern 
is very similar for other values of /3 and Rt- 



4.2 Density and pressure 

The density distribution associated with the flow is shown 
in Figure i.2. For a uniform atmosphere the density of the 



shocked gas decreases towards the cocoon to a value which is 
below that of the external density, similar to the behaviour 
of the gas behind a spherical bow shock. For /3 = 1 the 
density decreases slowly from the bow shock to the cocoon, 
again similar to the spherical case (Figure W), however, the 
density of the gas just behind the shock now decreases mono- 
tonically with increasing x because of the external density 
gradient. For /3 = 2 a similar trend with increasing x is seen, 
together with an increase in density from the bow shock to 
the cocoon. 

The pressure distribution of the gas between bow shock 
and contact discontinuity is shown in Figure 4.2, For a uni- 



form external density the pressure decreases monotonically 
from the hot spot towards the a;-axis, whereas for a decreas- 
ing density profile the pressure firstly decreases away from 
the hot spot then increases again towards the x-axis. For 
/3 = 2 the pressure close to the a:;-axis almost reaches the 
pressure in front of the hot spot. 

In the analysis presented in KA which led to the pre- 
diction of self-similar expansion of bow shock and cocoon 
we showed that it is a good approximation to assume a con- 
stant pressure distribution within the cocoon due to a high 
sound speed in this region. Across the contact surface pres- 
sure is of course continuous. The results shown in Figure 



4.2 indicate a pressure variation around the contact surface 
which is a minimum at some point between the hot spot 
and the a;-axis. In fact, within the cocoon there must be 
a backfiow (e.g. Norman et al. 1982) which is initially ac- 
celerated from the hot spot and must, given the assumed 
symmetry, be decelerated towards the symmetry plane of 
the a;-axis. For most values of /? this pressure variation is 
small and therefore we believe the expansion of the source 
will remain approximately self-similar and the results of KA 
and the analysis presented here will apply. For /3 close to 
2 the pressure variation becomes significant. The solution 
presented here will then be valid while the source has an ap- 
proximately elliptical shape, however in this case the source 
is unlikely to expand self-similarly. This is because one of the 
assumptions for self-similarity is the pressure equilibrium of 
the jets within the radio cocoon the pressure within which is 
assume to be uniform away from the hot spots. At least two 
possibilities for the evolution of such a source exist. Firstly, 



if the source retains approximate axial symmetry then as 
the high sound speed in the cocoon leads to pressure equal- 
ization, the excess pressure in the shocked gas will cause the 
cocoon to collapse resulting in a 'pinching off' of the cocoon; 
the bow shock will then not follow an elliptical coordinate 
surface, although subsequent expansion may still be approx- 
imately self-similar. Alternatively, the axial symmetry of the 
problem may not be retained leading to significant off-axis 
flow in the cocoon as is commonly observed. 

Note also that with the steeper rise of the pressure in 
direction of the core for large /3 the density is increasing as 
well in these regions. This implies that radiative cooling due 
to thermal bremsstrahlung becomes important and our as- 
sumption of adiabatic conditions may not be valid at this 
point. Any cooling will decrease the pressure and includ- 
ing radiative processes in the analysis should alleviate the 
problem of the steep pressure gradients towards the core. 
However, including these effects is beyond the scope of this 
paper and in the present form the existence of the pressure 
gradients described above means that the model is not fully 
self-consistent for large values of (3. 

A comparison of the density and pressure distributions 
for different aspect ratios of the bow sho ck b ut for the same 
external atmosphere is shown in Figure 4.2. The gradients 



in the density and the pressure distribution are greater for 
higher Rb- The variation in pressure is again seen to be 
smaller than in density. The stand-off distance between bow 
shock and cocoon at the hot spot. Ah, decreases for increas- 
ing Rb, while it stays almost constant for constant Rb but 
increasing /3 (see Table 0) . For large values of Rb the stand- 
off distance at the a:;-axis, A^,, decreases with increasing /3. 
We conclude that for /3 somewhat less than 2 self-similar 
expansion can occur with a flow distribution and source 
shape well approximated by the analysis presented here. For 
j3 close to 2 the source may pass through a phase in which 
the analysis of the present paper applies, but this does not 
represent a self-similar expansion and deviations from the 
assumptions of axial symmetry or elliptical bow shocks must 



4.3 Confinement of the cocoon 

In our analysis presented in KA the ratio of the pressure at 
the hot spot, Ph, and that of the cocoon, pc, determines the 
value of the dimensionless constant ci which describes the 
expansion of the cocoon. In the case of a cylindrical cocoon 
which is confined by the ram pressure of the receding IGM 
a simple estimate gives 
Ph 

Pc 

This estimate ignores the variation in external density 
along the outside of the cocoon and the small amount of 
backfiow. Taking pc to be the pressure at the contact dis- 
continuity half way between the core and the hot spot, we 
find that our results are well fitted by empirical relations of 
the form 



AR'' 



(16) 
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The assumption of ram pressure confinement of the cocoon, 
equation (hq) , therefore leads to an overestimate of the ratio 

Ph/Pc- 

Because of the self-similar expansion of the cocoons, the 
cocoon volume is given as csL^, where Lc is the length of 
the cocoon and C3 is a constant depending on the cocoon ge- 
ometry. The values of cs derived from the analysis presented 
here are given in Table ^ For the assumption of a cylindrical 
cocoon C3 equals 7r/(4 7?^); we find a maximal deviation of 
only 28% of this approximation from our results. 

The profile of the density distribution in the environ- 
ments of extragalactic radio sources is well fitted by a King 
(1972) profile. For distances from the centre of the density 
distribution greater than the core radius, Uo, the King profile 
is well fitted by the power law distribution, and for distances 
smaller than Oo a constant density po is a good approxima- 
tion. Consider a radio source which starts out in such a 
uniform density environment and expands beyond the core 
radius. KA argue that the ratio Ph/pc is determined by in- 
trinsic properties of the jet since the cocoon and the jet are 
assumed to be in pressure balance. In this case equations 
( [l7[ ) imply that the cocoon of such a source will become nar- 
rower when it emerges from the uniform density regime into 
the region with a decreasing density gradient. This result 
is in agreement with the observation of Black (1992) who 
notes that in his sample shorter sources tend to be 'fatter'. 



5 X-RAY EMISSION FROM THE HEATED 
IGM 

Assuming the shocked IGM between bow shock and contact 
discontinuity is an ideal gas consisting entirely of ionised 
hydrogen the absolute temperature at any given point in 
the flow is 






(18) 



where irip is the mass of a proton and k is the Boltzmann 
constant. Note, that the temperature of the gas in the flow 
does not depend on the temperature of the environment the 
source is expanding into. 

Figure n shows a plot of the ratio P/R for various at- 
mospheres. The difi'erences between the temperature distri- 
bution for different /3 are relatively small. The thin shell of 
rarefied gas around the cocoon in an uniform atmosphere 
is slightly hotter than its surroundings and the overall tem- 
perature of the gas in the fiow decreases somewhat for in- 
creasing /3. A comparison of models for various aspect ratios 
(Figure 0, left column) shows an increasing prominence of 
the hot spot region with increasing Rb which can be ex- 
plained by the faster expansion speed of the bow shock at 
the hot spot with respect to the expansion perpendicular to 
the jet axis. Again the temperature scale is very similar in 
all four panels. 

The specific emissivity of ionised hydrogen at tempera- 
ture T due to thermal Bremsstrahlung is (e.g. Shu, 1991) 



= 7 X 10 



-51 rie Up 



Vt 



w 

i3Hz 



where n^ and n^ are the number densities of the electrons 
and of the protons respectively. When observed with a finite 
bandwidth, ui -^ U2, the integrated X-ray emissivity is given 
by 



e = 5.6 X 10" ^' °> L^^- v^Wp 



ai 

-4.0x10"'' a\ L^ '^ V ^ 



m3 



with 



a\ 



5-/? 



/3\3 



Ci \Po0.o) ^ Qo 



(20) 



(21) 



where we have used equation (bI). To obtain the X-ray sur- 
face brightness this expression has to be integrated along the 
line of sight through the flow. The dependence of e on the ad- 
vance speed of the bow shock. Lb, is non-linear. To proceed 
we therefore have to choose specific values for the external 
density profile, the jet power, the linear size of the cocoon, 
Lj and the observing frequency. Assuming po = 5 x lO"'^^ 



kg/m'^ 



1 kpc, typical values for the gas surround- 



(19) 



ing isolated galaxies (Canizares ei al. 1987), Qo ~ 10"^^ W 
(Rawlings & Saunders 1991), Lj = 100 kpc and an observ- 
ing band from ui =0.1 keV to !/2 = 2 keV (appropriate for 
the High Resolution Imager, HRI , of the ROSAT satellite) 
we obtain the results shown in Figure 0. Here we have as- 
sumed that the source is rotationally symmetric about the 
2-axis, that the source lies exactly in the plane of the sky and 
we have neglected the cosmological redshift of the observing 
frequency. 

The appearance of the fiow in X-ray emission varies 
significantly with j3. For a uniform density environment the 
hot spot region dominates the emission, while for higher val- 
ues of P the region close to the i-axis is brightest and the 
emission is concentrated towards the edge of the cocoon. 
Together with the observation made earlier that the tem- 
perature variation is not very large within the flow it is now 
evident that the X-ray surface brightness for the frequency 
band chosen is a good tracer of the density distribution. X- 
ray observations of FRII sources with high resolution can 
therefore constrain the shape of the density distribution of 
the atmosphere the source is expanding into. A similar result 
was found by Clarke et al. (1997) from numerical simulations 
of a cocoon-bow shock structure around a jet. For smaller 
aspect ratios the X-ray surface brightness is more diffuse 
(see Figure pi right column) , but the general features of the 
appearance of the source remain the same. Note however the 
difference in the ages of the sources of the same linear size 
but of different geometry and/or environment. 

Heinz et al. (1998) introduce observational diagnostics 
for X-ray observations of extragalactic radio sources. These 
are based on their assumption of spherical geometry for 
the bow shock and cocoon which implies a uniform surface 
brightness at all points within the flow at equal distance 
from the centre of the source. From the analysis presented 
here it is clear that the X-ray surface brightness of the flow 
between bow shock and cocoon can vary considerably within 
the fiow region if the source geometry is elongated and/or 
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the density gradient of the external medium is steep. In these 
cases it is not straightforward to infer the absolute values of 
source or environment parameters from X-ray observations 
alone. 

The total radiation due to thermal bremsstrahlung, 
dE/dt, can be found by setting !/i = and i^2 — > oo in 
equation (GOl) and integrating the resulting emissivity over 
the volume of the flow. For almost all the cases investigated 
here we find the total luminosity to be significantly less than 
10% of the jet power. We therefore conclude that our adia- 
batic treatment of the flow between bow shock and cocoon 
is justified. 

A rough estimate for the local cooling time of the 
shocked IGM can be obtained by dividing the local energy 
density, -P(T?,f), by dE/dt. For all cases presented in this 
paper we find cooling times of the order of 10% to a few 
thousand times the Hubble time which is much greater than 
the inferred life times of radio sources. The presence of a 
radio source therefore infiuences the evolution of its envi- 
ronment far beyond the limited life time of the source itself. 



5.1 Application of the model to Cygnus A and 3C 
356 

The high resolution ROSAT map of CariUi ei al. (1994) has 
sufficient spatial resolution to allow a comparison of source 
features observed in the radio at low frequency with those 
observed in X-rays. For the density distribution of the mate- 
rial surrounding Cygnus A we find po = 1.4 x lO"'^^ kg/m'', 
ao — 25.6 kpc and l3 = 1.4 from the fitting of the X-ray 
surface brightness distribution expected from a King profile 
(n = 1, Ho = 50 km s~^ Mpc~^) to the data presented by 
CariUi et al. (1994). Note, that this is different from the fit 
obtained by these authors. Arnaud et al. (1987) determine 
the temperature of the gas to be 4 x 10^ K. CariUi et al. 
(1994) also present a 347 MHz VLA map of Cygnus A from 
which it is possible to determine the linear size and the as- 
pect ratio of the cocoon. Assuming that Cygnus A is lying 
exactly in the plane of the sky we find Lj « 105 kpc and 
R ~ 2; the difference in the lengths of the two cocoons is 
roughly 15%. The results of an integration of the hydrody- 
namical equations describing the flow between bow shock 
and cocoon for these parameters is given in Table m The 
aspect ratio of the bow shock of Cygnus A implied by this 
calculation is 1.6. 

Using the model of Kaiser, Dennett-Thorpe & Alexan- 
der (1997) for the radio emission of FRII sources we have 
calculated the radio spectrum of Cygnus A for the given 
parameters. The model is in good agreement with the ob- 
served radio luminosities in the Gigahertz range presented 
by Baars et al. (1977) for a jet power of Qo = 2 x 10^" W. 
Using the value of Ph/Pc = 4.1 and the model of KA we find 
an age of 2.9 x 10^ years for Cygnus A. The present hot spot 
advance speed implied by this is 9.7 x 10~^ c. Both results 
agree well with those of CariUi et al. (1991). 

After subtracting the expected X-ray emission due to 
thermal bremsstrahlung of a King profile with the parame- 
ters quoted from the map presented by CariUi et al. (1994), 
the excess or deficit of the emission due to the presence of the 



radio source becomes visible. The most interesting feature 
for this analysis are the emission enhancements to the sides 
and just outside the cocoons because they must be caused 
by the shocked IGM between bow shock and cocoon. For the 
eastern cocoon the peak in the residual emission lies towards 
the core of the radio structure and the excess measured af- 
ter the King profile has been subtracted is 0.21 counts per 
second per 0.5 arcsecond pixel with the HRI in the band 
between 0.1 keV and 2 keV. This implies an X-ray emission 
of the gas in this region corresponding to 0.27 counts per 
second per pixel with the HRI because the measured count 
rate along a line of sight through the cocoon region, which 
we assume to be free of X-ray emitting gas, is 0.06 counts 
per second per pixel less than that expected from the King 
profile. For the western cocoon the situation is less clear and 
the emission enhancement may not be significant. 

The surface brightness map for the case of Cygnus A is 
very similar to the one shown on the bottom right of Figure 
pi Qualitatively the predictions of the numerical calculation 
presented here agree well with the observations. The emis- 
sion enhancements are at the correct positions, close to the 
core of the source just outside the cocoon. However, a closer 
examination of the temperature distribution reveals some se- 
vere problems with the application of the results found in the 
previous section to the case of Cygnus A. The value of P/R 
at the point of highest surface brightness is roughly 0.012. 
Together with the derived advance speed of the hot spot we 
find a temperature of 1.9 x lO'' K at this point. This is well 
below the temperature of the King profile assumed for the 
gas which surrounds the bow shock. Indeed, even the model 
temperature just behind the bow shock close to the core of 
the source is below the temperature of the unshocked gas. 
This is caused by the assumption of strong shock conditions 
along the entire bow shock which may not be applicable to 
the case of Cygnus A. 

Using the temperature of the ambient gas the peak 
in the X-ray surface brightness predicted by our model is 
4 X 10~^^ W m~^ per 0.5 arcsecond pixel which corresponds 
to 0.012 counts per second per pixel with the HRI. For the 
conversion we have used the Internet version of PIMMSH 
with a column density for the galactic medium of 3 x 10 
m~^ towards Cygnus A (Clarke et al. 1997). This discrep- 
ancy may be explained by the existence of cold, dense clumps 
of gas embedded in the otherwise smooth distribution of gas 
in galaxy clusters (Fabian et al. 1994). The X-ray emission 
of a volume of gas after the bow shock has passed through 
it is proportional to p^My/T, where p is the pre-shock den- 
sity of the gas, M the Mach number of the shock within 
the gas and T the temperature. For ideal gas conditions 
the Mach number is proportional to v/VT, where v is the 
velocity of the shock with respect to the gas. The ratio of 
the X-ray luminosity of a cold, dense cloud and that of hot 
IGM with the same volume is then given by PcVc/{pxVx) 
, where subscripts c and x are used for the cloud material 
and the hot IGM respectively. The cold clouds are small and 



■f PIMMS was programmed by K. Mukai at the High Energy 
A strophysics Science Archive Research Center of NAS A, available 
at http://heasarc.gsfc. nasa.gov/Tools/w3pimms.htm\ 
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will therefore quickly re-establish pressure equilibrium with 
the shocked hot IGM surrounding it after the bow shock 
has passed through it (e.g. McKee & Cowie 1975). We can 
therefore set the veloci ty of the bow shock within the cloud 
material to ~ \fpxlpcVx ■ With this we find that the ratio of 
the X-ray luminosity of a cold cloud and that of the hot IGM 
occupying the same volume is given by (pc/px) ■ The peak 
in the X-ray emission observed in the eastern lobe of Cygnus 
A can thus be explained if the gas in this region is about 5 
times denser than the otherwise smoothly distributed IGM. 

The striking symmetry of the two emission peaks on 
either side of the base of the eastern cocoon is difficult to 
explain in this scenario. Note however, that there is another 
peak in the X-ray emission situated to the south-west of the 
core well inside the western cocoon. This could be another 
clump of cold material with similar properties to the two dis- 
cussed above, which is observed in projection and therefore 
appears to lie within the cocoon. 

Crawford & Fabian (1996) report the detection of ex- 
tended X-ray emission in the vicinity of the FRII radio 
galaxy 3C 356. They argue that the emission is the signature 
of an extended IGM in a cluster around 3C 356 which makes 
this radio galaxy the most distant object (z=1.079) in which 
cluster gas has been detected directly. The resolution of the 
X-ray map of this object is insufficient for any quantitative 
study. However, the observed extension is aligned with the 
northern cocoon of 3C356 which suggests that the X-ray 
emission of the IGM surrounding the northern cocoon may 
be enhanced by the compression and heating of this material 
by the bow shock of the radio source, similar to the case of 
Cygnus A. 



5.2 An empirical model for the X-ray emission 

It is straightforward to calculate the total X-ray luminosity 
of the shocked IGM between the bow shock and the cocoon 
for a given frequency and bandwidth, Lx, predicted by the 
model for a given set of source parameters by integrating 
over the surface area of the bow shock. However, this re- 
quires the numerical calculation to be performed for this 
particular set of parameters. In the following we develop 
an analytical approximation for the X-ray luminosity using 
empirical relations of source parameters derived from the 
numerical calculations. 

The total mass of the external gas swept up by the bow 
shock during the life time of the radio source, M, is given 
by 



M = 



where 



Po 



r/2 



dVt 



■ 4 7r/9offlfib 



'"' . , . ^ sinh^ ri + sin^ ^ ,^ , 

smh r? sm ^ — . d(, d-q. 

( cos^ ^ -I- smh 77J 



(22) 



(23) 



If we assume this mass to be uniformly distributed 
over the volume occupied by the flow, the density in this 
region is found to be p = M/{Vb ~ 14) = M/{a2Ll), 
where Vt and Vc are the volume enclosed by the bow 



shock and that of the cocoon respectively. By definition 
Vb — 2tvLI (Co sinh^ rjb — l/3tanh^ rytj and Vc = C3L5. 

We define T as the average of the temperature of the 
gas in the flow and since the temperature variations within 
the flow material are small, local deviations from T will be 
small as well. For an ideal gas T = rup/k LI P/R. 

With the two assumptions of uniform density and uni- 
form temperature in the flow we find from equation (M) 



5.2 X 10" p- 



9.1 X 10' 



/r {Vb - Vc] 

[po fflf ) 



^/T 



■r-^li 



a2 

T-7 „2 



4-2^j 



P/R 



l/P/R 



(24) 



This expression depends on ci, P/R and I/a2 which are 
all functions of the exponent of the density profile of the 
unshocked IGM, /3, and the geometry of the FRII source 
defined by the aspect ratios R or Rb. Empirical expressions 
for these functions can be derived from the results of the 
numerical calculation. We find 





' 0.17i?-'-^" 
< 0.13 i?-^-^* 
. 0.095 i?-^-^^ 
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0.75 7?°-5" 
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1.35 Rt"^^ ■,f3 = 
1.18 R^b""^ ;/? = 
0.85 Rl^^ -,13 = 



1 
2 





(25) 



Figure p^ shows a comparison of equation iVM using 
the empirical relations above, with the results of the numer- 
ical calculation for po = 5 x 10~^^ kg m~"^, Oo = 1 kpc, 
Qo = 10=^^ W, 1^1 = 0.1 kev and 1/2 = 2 keV. For /3 > the 
empirical model agrees well with the numerical results. The 
deviation for large linear sizes and /3 = is acceptable in 
many cases since it is unlikely that a uniform density region 
extends to very large distance from the centre of any object 
hosting a radio source. Note however, that the fit of the em- 
pirical model worsens for higher observing frequencies since 
the exponential part of equation (pOl), which is not well fitted 
by the empirical formula, then dominates this expression. 



6 CONCLUSIONS 

We presented a numerical integration of the hydrodynamical 
equations governing the self-similar gas flow behind a strong 
shock in two dimensions assuming the shock to have an ellip- 
soidal geometry. The density distribution in the unshocked 
environment of the shock is modeled with a power law. Ap- 
plying this model to the bow shocks surrounding the cocoons 
of extragalactic radio sources we found that there should be 
a pressure gradient within the cocoons of sources located in 
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non-uniform environments which is consistent with the ex- 
istence of backflow within the cocoon. In extreme density 
profiles non-axiaUy symmetric flow or pinching-off of the co- 
coon is very likely to occur. Furthermore, the expansion of 
the cocoon and bow shock may not be self-similar. In these 
cases (/3 ~ 2) the model is not self-consistent. We also find 
that the assumption of ram pressure confinement of the co- 
coon perpendicular to the jet axis leads to an overestimate 
of the ratio of the pressure in front of the hot spot and the 
pressure in the cocoon. 

From the properties of the gas in the flow we calculated 
the X-ray surface brightness predicted by the model. Be- 
cause of the small variations of the temperature in the flow, 
the surface brightness was found to be a good tracer of the 
gas density. The appearance of radio sources in X-rays is 
found to vary significantly with the properties of the exter- 
nal medium. Simple diagnostic tools have to be refined to 
account for the elongation of the bow shock of most sources. 
The cooling times of the IGM shocked by the bow shock are 
much longer than the expected life time of the radio sources 
causing the bow shock. The evolution of any concentration 
of matter in the universe will therefore be significantly infiu- 
enced by the presence of a radio source in its centre. This is 
the case even if the activity time scale of the source is much 
shorter than the evolution time scale of its surroundings. 

We compared the predictions of our model with the X- 
ray map of Cygnus A of Carilli et al. (1994). Although our 
model predicts the strongest X-ray emission at the position 
where it is observed, on the outside of the cocoon closest 
to the core of the source, the model can not reproduce the 
observed luminosity. The emission may be caused by the 
presence of overdense clumps of gas embedded in the other- 
wise smooth IGM. 

The enhanced X-ray emission of the shocked IGM be- 
hind the bow shock of a strong radio source may also explain 
the observed extended X-ray emission around 3C 356. 

Based on the numerical calculation we also presented 
an empirical model for the total X-ray emission expected 
from the shocked IGM surrounding the cocoons of powerful 
FRII sources. 

Despite the inconsistencies within the model that arise 
for large values of (5 we believe that this analysis is an impor- 
tant first step in determining the influence of powerful radio 
sources on the properties and evolution of their gaseous en- 
vironments. 
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APPENDIX 

A.l Comoving prolate spheroidal coordinates 

A. 1.1 Fixed coordinates 

The Cartesian coordinates of any point in space can be re- 
lated to their prolate spheroidal counterparts via the expres- 



z — c cosh 77 cos^, 

X — c sinh 77 sin ^ sin (j>, 

y — c sinh 77 sin ^ cos (j>, 



(26) 



where c is a constant describin g the e ccentricity of the ellip- 
soids of constant 77 (see Figure A. 1.1 ). Since we assume that 
the bow shock is of ellipsoidal shape with the z-axis, the jet 
axis, being the axis of rotational symmetry the problem will 
be independent of (j) and we seek a solution (without loss of 
generality) in the xz-plane, i.e. (p = 7r/2. The unit vectors 
with respect to rj and ^ are 
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e,, = — (sinh 77 cos^ e^ + cosh 77 sin^ e^,) , 

Cr 

e^ = — (— cosh7; sin^ e^ + sinh77 cos^ e^,) , (27) 
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with 



Cr — v sinh'^ ri + sin'^ ^ 

It immediately follows that 



(28) 



Grj — Cg G^ , Gfj — Ce 6^ 



9 . _ _. 
where 



9 . 



9?? 



ej — Cg e,, 



(29) 



sinhr; cosh 77 



sin^ cos^ 



(30) 



Using equations (E9) the r] and ^ components of the 
V-operator in this coordinate system can be found. 



V„ = 



1 9 _ 1 
ccr dr] cc, 



V' 



1 9 



Ve = — — EE — V^ 

CCr Ot, CCr 



(31) 



The divergence of a vector, V • V, in prolate spheroidal co- 
ordinates assuming that V does not depend on (j) can be 
written in the form 



CCrV-V = — V„ + (ce+coth?7)V^ 



(32) 



A. 1.2 Comoving coordinates 

By setting c = CoLt it is possible to transform the station- 
ary coordinate system discussed in the previous section into 
comoving coordinates. From equations (G6|) it is clear that 
if c is a function of time, rj and ^ must also depend on time. 
Taking the derivative of (feq) with respect to t one finds 






(33) 



where a dot denotes a time derivative. This also means that 
all partial derivatives with respect to time in equations (W) 
are replaced according to 

d_ ^ d_ 
dt ^ dt 



drj 95 



(34) 



Note, that the partial derivatives with respect to the spatial 
coordinates are unchanged and we therefore continue to use 
the symbols, rj and ^, for these spatial coordinates. 
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Table 1. The results of the numerical integration. fS is the exponent of the power law for the external density distribution, ij(, is the 
aspect ratio of the bow shock, R is the aspect ratio of the cocoon, Ph/Pc is the ratio of the pressure at the hot spot and the pressure in 
the cocoon defined as the pressure at the contact discontinuity half way between the core of the radio source and the tip of the cocoon, 
?7f, is the elliptical coordinate surface coinciding with the bow shock, »?^ is the elliptical coordinate at which the flow reaches the cocoon 
in front of the hot spot, rjc is the elliptical coordinate at which the flow reaches the point on the contact discontinuity half way between 
the core and the tip of the cocoon, rj^ is the elliptical coordinate at which the flow reaches the cocoon on the x-axis, A^ is the stand-off 
distance of the bow shock from the contact discontinuity at the hot spot, A^; is the stand-off distance of the bow shock from the cocoon 
at the X-axis and C3 is the volume constant. 

Rb R Ph/Pc Vb Vh Vc Vx Aft Ax C3 







1.0 


1.221 


3.518 


0.6585 


0.5775 


0.5386 


0.5373 


0.0435 


0.1160 


0.4153 






2.0 


2.679 


16.92 


0.2971 


0.2581 


0.2217 


0.2217 


0.0105 


0.0741 


0.0940 


/3 = 





3.0 


4.125 


39.80 


0.1949 


0.1699 


0.1415 


0.1415 


0.0045 


0.0528 


0.0403 






4.0 


5.559 


71.63 


0.1454 


0.1268 


0.1044 


0.1044 


0.0025 


0.0406 


0.0223 






5.0 


6.985 


111.8 


0.1160 


0.1022 


0.0828 


0.0828 


0.0015 


0.0329 


0.0146 






1.0 


1.217 


2.384 


0.6585 


0.5715 


0.5386 


0.5360 


0.0465 


0.1172 


0.4153 






2.0 


2.579 


7.973 


0.2971 


0.2581 


0.2300 


0.2276 


0.0105 


0.0683 


0.1015 


/3 = 


: 1 


3.0 


3.858 


16.67 


0.1949 


0.1699 


0.1512 


0.1477 


0.0045 


0.0466 


0.0461 






4.0 


4.997 


28.34 


0.1454 


0.1268 


0.1160 


0.1142 


0.0025 


0.0308 


0.0276 






5.0 


6.297 


43.34 


0.1160 


0.1022 


0.0919 


0.0919 


0.0015 


0.0239 


0.0174 






1.0 


1.211 


1.509 


0.6585 
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0.5400 


0.5321 


0.0485 


0.1209 
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0.0695 
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0.0045 


0.0470 


0.0527 






4.0 


4.758 


11.96 


0.1454 


0.1268 


0.1218 


0.1218 


0.0025 


0.0232 


0.0305 






5.0 


5.895 


17.99 


0.1160 


0.1022 


0.0981 


0.0981 


0.0015 


0.0176 


0.0199 














Cygnus 
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/3 = 


: 1.4 


1.6 


2.031 


4.099 


0.3779 


0.3256 


0.2970 


0.2917 


0.0175 


0.0842 


0.1610 
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Figure 1: The flow behind a spherical bow shock in various atmospheres. The top panel shows the velocity of the flow, 
the bottom panels the density and the pressure. Solid curves; l3 — 0, dashed curves; l3 — 1 and dotted curves: /3 = 2. 



Figure 2: The velocity field of the flow between bow shock and contact discontinuity in the frame of the comoving 
coordinate system for /3 = and Rb — 1. The bow shock deflnes the upper edge, the contact discontinuity the lower edge of 
the field plotted. The AGN which drives the jet is located in this plot at the point {z = 0, x = 0) and the tip of the bow 
shock is at (2 = 0, a; = 1). 



Figure 3: Details of the flow pattern shown in Figure 4.1. Left; the flow close to the hot spot. Right; the flow near the 



a;-axis. Note, that the vectors representing the flow velocity have the same length scale in both plots. 



Figure 4: Density profile of the flow between bow shock and contact discontinuity. Top: Uniform atmosphere, (3 = 0. 
Bottom left: P — 1. Bottom right: /3 = 2. For all plots Rt — 3. The gray scales are normalised such that the density just ahead 
of the tip of the bow shock is the same in each panel. 



Figure 5: Pressure profile of the fiow between bow shock and contact discontinuity. Top; Uniform atmosphere, /3 = 0. 
Bottom left: /3 = 1. Bottom right; /3 = 2. For aU plots Rb = 3. 



Figure 6: Comparison of the density and pressure distribution for various aspect ratios of the bow shock. Left column: 
density, right column; pressure. From top to bottom; Rb — 1, Rb = 2, Rb = 4 and Rb — 5. For all plots /3 = 1. 



Figure 7: Temperature profile of the fiow between bow shock and contact discontinuity. Top: Uniform atmosphere, /3 = 0. 
Bottom left; (3 = 1. Bottom right; f3 = 2. For aU plots Rb = 3. 



Figure 8: Comparison of the temperature and X-ray surface brightness distribution for various aspect ratios of the bow 
shock. Left column: temperature, right column: surface brightness. From top to bottom: Rb = 1, Rb = 2, Rb = 4 and Rb = 5. 
For all plots f3 = 1. For the chosen parameters (see text) the ages of the radio sources in the right column are 3.0 x 10^ years, 
1.5 X 10^ years, 8.8 x lO" years and 7.4 x 10^ years respectively from top to bottom. 
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Heating of the IGM by FRII radio sources 13 

Figure 9: Surface brightness profile of the fiow between bow shock and contact discontinuity. Top: Uniform atmosphere, 
/3 = 0. Bottom left: /3 = 1. Bottom right: P — 2. For all plots Rb = 3. See text for the parameters of the external atmosphere 
and the jet power. The resulting ages of the radio sources are: Top: 3.7 x 10^ years, bottom left: 1.1 x 10^ years and bottom 
right: 3.5 x 10^ years. 

Figure 10: Comparison of the empirical model for the X-ray emission with the numerical calculation. Individual symbols 
are results from the numerical calculation while curves represent the empirical model. Top panel: Rt — 1.0, bottom left: 
Rb = 3 and bottom right: Rt = 5. Solid lines and circles: /3 = 0, dashed lines and crosses: /3 = 1 and dotted lines and squares: 
(3 = 2. See text for other model parameter. 

Figure 11: Prolate spheroidal coordinates for c = 0.6 and rotation about the z-axis. 
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